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Abstract 


Background 


A highly pathogenic human coronavirus (CoV), Middle East respiratory syndrome coronavirus 
(MERS-CoV), has emerged in Jeddah and other places in Saudi Arabia, and has quickly 
spread to European and Asian countries since September 2012. Up to the 1$ October 2015 it 
has infected at least 1593 people with a global fatality rate of about 35%. Studies to understand 
the virus are necessary and urgent. In the present study, MERS-CoV main protease (MP"°) is 
expressed; the dimerization of the protein and its relationship to catalysis are investigated. 


Methods and Results 


The crystal structure of MERS-CoV MP" indicates that it shares a similar scaffold to that of 
other coronaviral MP and consists of chymotrypsin-like domains | and II and a helical 
domain Ill of five helices. Analytical ultracentrifugation analysis demonstrated that MERS- 
CoV MP"? undergoes a monomer to dimer conversion in the presence of a peptide substrate. 
Glu169 is a key residue and plays a dual role in both dimerization and catalysis. The muta- 
genesis of other residues found on the dimerization interface indicate that dimerization of 
MERS-CoV MP" is required for its catalytic activity. One mutation, M298R, resulted in a sta- 
ble dimer with a higher level of proteolytic activity than the wild-type enzyme. 


Conclusions 


MERS-CoV MP" shows substrate-induced dimerization and potent proteolytic activity. A 
critical assessment of the residues important to these processes provides insights into the 
correlation between dimerization and catalysis within the coronaviral MP? family. 


Introduction 


A highly pathogenic human coronavirus (CoV)’, Middle East respiratory syndrome coronavi- 
rus (MERS-CoV), emerged in Jeddah and other places in Saudi Arabia in September 2012 and 
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spread to some European, African and Asian countries in recent years [1-3]. The virus causes 
symptoms similar to severe acute respiratory syndrome coronavirus (SARS-CoV), but an addi- 
tional component is also involved, namely acute renal failure [4]. Up to the 1** October 2015, 
1593 people had been infected with MERS and this has led to 568 reported deaths globally 
(World Health Organization, global alert and response, http://www.who.int/csr/don/ 
01-october-2015-mers-jordan/en/). A serological study of major livestock in Saudi Arabia sug- 
gested dromedary camels may be an original reservoir [5, 6]; although recent studies have 
identified bats may also as the suspected host as a number of bat CoVs show high sequence 
similarity to SARS-CoV and MERS-CoV [7, 8] and a bat-CoV was discovered that readily 
infect human cells using ACE2 as the receptor [9]. Nevertheless, human-to-human transmis- 
sion of MERS-CoV has been confirmed. This May, an infected traveler from the Middle East 
region returned to his home country of the Republic of Korea and caused an outbreak that 
revolved around health facilities (http://www.who.int/csr/don/07-july-2015-mers-korea/en/). 
These findings indicate that it is possible for the virus to spread globally and as such it poses a 
significant threat to world health and the world economy in general. Therefore studies that 
help to understand this virus and aid the development of antiviral drugs or other therapies are 
important. 

Similar to other CoVs, after entering the host cells, the nonstructural polyproteins (ppla 
and pplab) of MERS-CoV are synthesized and then cleaved by two coronaviral proteases, a 
main protease (M?"®) (EC 3.4.22.69) and a papain-like protease (EC 3.4.22.46) [10]. This cleav- 
age is considered to be a leading process which is required for viral maturation [11-14]. The 
MERS-CoV MP"®, namely nsp5 of the pp1a proteins (residue 3248-3553), has been identified 
[15]. Like other M?"°, there is a catalytic dyad that consists of a His residue and a Cys residue 
[16-20]. Sequence alignment suggests that MERS-CoV MP", in a similar manner to other 
known MP", has a chymotrypsin-like architecture consisting of a catalytic core (domain I and 
II) and a helical domain III; sequence identities of CoV M?"° protein range from 50% to 80% 
(S1 Fig). Recently the crystal structure of an inactive MERS-CoV MP" C148A mutant has con- 
firmed this similarity and the results also suggest that this protein forms a dimer [15]. Further- 
more, based on the identification of eleven canonical cleavage sites, the MERS-CoV MP"° 
should be able to recognize and cleave at the (L/M)-Q-|-(A/S) conserved sequence, which is 
essential for most CoV M?"°-mediated processing [21, 22]. However, up to now, the correlation 
between the protein’s structure and the catalytic process remains unclear. 

In the present study we expressed and purified the MERS-CoV MP"® using the authentic N- 
terminus via an Escherichia coli system. The crystal structure of the MERS-CoV MP?” at 3.0-A 
resolution is reported. The quaternary structural changes in the MERS-CoV MP" in the 
absence and presence of peptide substrates were investigated by analytical ultracentrifugation 
(AUC). The results of kinetic activity assays indicated that MERS-CoV MP?” exhibits potent 
proteolytic activity that is associated with a pattern of cooperativity. Some critical residues for 
dimerization of the protein and catalysis by the protein were verified by site-directed mutagen- 
esis. The present studies provide a foundation for an understanding of the mechanism that 
controls the monomer-dimer switch at work in MERS-CoV MP”. 


Materials and Methods 
Expression Plasmid Construction 


The sequence of the MERS-CoV MP"° (GenBank accession number AHC74086; polyprotein 
residues 3248-3553) was synthesized (MDBio Inc.), digested by NdeI-XholI and then inserted 
into the vector pET-28a(+) (Novagen). In this construct, the 6 x His tag is retained at the 
N-terminus. To remove the fusion tag and generate an authentic N-terminus for protein 
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purification, the codons of the thrombin cutting recognition sequence and a Ndel cutting site 
were removed and then inserted the codons of Leu-Arg-Leu-Lys-Gly-Gly into the above vector. 
The forward primer sequence for site-directed mutagenesis was 5’--CATCACAGCAGCGGCCT 
GCGTCTGAAAGGCGGCAGCGGTTTGGTGAAAATG-3’ and the reverse primer was 5’-CATTTT 
CACCAAACCGCTGCCGCCTTTC AGACGCAGGCCGCTGCTGTGATG-3’. The reading frame of 
the final plasmid was confirmed by sequencing. 


Expression and Purification of MERS-CoV MP? 


The expression vector was transformed into E. coli BL21 (DE3) cells (Novagen). Cultures were 
grown in 0.8 liters of LB medium at 37°C for 4 h, induced with 0.4 mM isopropyl-B-p-thioga- 
lactopyranoside, and then incubated overnight at 20°C. After centrifuging at 6,000 x g at 4°C 
for 15 min, the cell pellets were resuspended in lysis buffer (20 mM Tris, pH 8.5, 250 mM 
NaCl, 5% glycerol, 0.2% Triton X-100, and 2 mM B-mercaptoethanol) and then lysed by soni- 
cation. The crude extract was then centrifuged at 12,000 x g at 4°C for 25 min to remove the 
insoluble pellet. Next the supernatant was incubated with 1-ml Ni-NTA beads at 4°C for 1h 
and then loaded onto an empty column. After allowing the supernatant to flow through, the 
beads were washed with washing buffer (20 mM Tris, pH 8.5, 250 mM NaCl, 8 mM imidazole, 
and 2 mM £-mercaptoethanol). The SARS-CoV papain-like protease [12] (1 mg in 100 mM 
phosphate buffer (pH 6.5)) was then added and incubated for 3 h. The SARS-CoV papain-like 
protease digestion, which removed the 6 x His tag and Leu-Arg-Leu-Lys-Gly-Gly fragment, 
resulted in a native protein product with an authentic N-terminus. The digest was allowed to 
flow through and then loaded onto a S-100 gel-filtration column (GE Healthcare) equilibrated 
with running buffer (20 mM Tris, pH 8.5, 100 mM NaCl, and 2 mM dithiothreitol). The purity 
of the fractions collected was analyzed by SDS-PAGE and the protein was concentrated to 30 
mg/ml by Amicon Ultra-4 10-kDa centrifugal filter (Millipore). 


Protein Crystallography 


Crystals of the MERS-CoV MP"° were obtained at 295 K by the sitting-drop vapor-diffusion 
method. The protein solution was set up at 5 mg/ml and the reservoir solution consisted of 0.1 
M Tris, pH 8.4, 15% (w/v) PEG 4000 and 0.2 M sodium acetate. Clusters of needle crystals 
appeared in 2 days and were used for micro-seeding. Single cystals of rectangle shape and with 
dimensions of 0.3-0.5 mm were obtained in less than a week. All crystals were cryoprotected in 
the reservoir solution with 15% glycerol and were flash-cooled in liquid nitrogen. 


Data collection, structure determination and refinement 


X-ray diffraction data were collected at 100 K on the SPXF beamline 13C1 at the National Syn- 
chrotron Radiation Research Center, Taiwan, ROC, using a ADSC Quantum-315r CCD detec- 
tor (X-ray wavelength of 0.976 A). The diffraction images were processed and scaled using the 
HKL-2000 package [23]. The structure was solved by the molecular replacement method by 
Phaser [24] using the structure of SARS-CoV MP"° R298A mutant (PDB entry 4hi3; [25]) as 
the search model. Manual rebuilding of the structure model was performed using Coot [26]. 
Structure refinement was carried out using REFMAC [27]. The data-processing and refinement 
has been deposited in the Protein Data Bank (PDB entry 5c3n). 


Steady-State Kinetic Analysis 


The colorimetry-based peptide substrate, TSA VLQ-para-nitroanilide (TQ6-pNA) (purity 95- 
99% by HPLC; GL Biochem Ltd, Shanghai, China), was used to measure the proteolytic activity 


PLOS ONE | DOI:10.1371/journal.pone.0144865 December 14, 2015 3/18 


®-PLOS | ove 


Dimerization and Catalysis of MERS-CoV Main Protease 


of MERS-CoV M?"® and its mutants throughout the course of the study as described previously 
[25, 28]. This substrate is cleaved at the GlIn-pNA bond to release free pNA, resulting in an 
increase in absorbance at 405 nm. The absorbance at 405 nm was continuously monitored 
using a Jasco V-550 UV/VIS spectrophotometer. The protease activity assay was performed in 
10 mM phosphate (pH 7.6) at 30°C. The substrate stock solution was 1600 uM and the working 
concentrations were from 25 to 1200 uM. In the substrate titration assay, the concentration of 
MERS-CoV MP*° and its mutants, V4R, T126S, E169A, M298R and T126S/M298R was 0.3, 
0.4, 0.7, 1.2, 0.15 and 0.26 uM, respectively, while that of SARS-CoV MP"° was 1.1 uM. Steady 
state enzyme kinetic parameters were obtained by fitting the initial velocity (vp) data to the 
Michaelis-Menten Eq (1) 


Yo >a 3. fer (1) 


where k.. is the catalytic constant, [E] is the enzyme concentration, [S] is the substrate concen- 
tration and K,, is the Michaelis constant of the substrate. The program SigmaPlot (Systat Soft- 
ware, Inc., Richmond, CA) was used for the data analysis. 

To assess the cooperativity effect, the kinetic parameters were obtained by fitting the initial 
velocities to the Hill Eq (2) 


— ku lEI(S]" 


~ K+ [sy] @) 


Vo 


where K’ is a constant that is related to the dissociation constant and h is the Hill constant. 


Analytical ultracentrifugation analysis 


AUC was performed on a XL-A analytical ultracentrifuge (Beckman Coulter) using an An-50 
Ti rotor [11, 12, 25, 28-30]. The sedimentation velocity experiments were carried out using a 
double-sector epon charcoal-filled centerpiece at 20°C with a rotor speed of 42,000 rpm. Pro- 
tein solutions of 0.05 to 0.5 mg/ml (330 ul) and reference (370 ul) solutions, both containing 
D,O, were loaded into the centerpiece. The absorbance at 280 nm was monitored in a continu- 
ous mode with a time interval of 300 s and a step size of 0.003 cm. Multiple scans at different 
time intervals were then fitted to a continuous c(s) distribution model using the SEDFIT pro- 
gram [31]. Additionally, the results with the various different protein concentrations were 
globally fitted to a monomer-dimer self-association model using the SEDPHAT program to 
calculate the dissociation constant (Kg) [32]. 

To measure the substrate-induced dimerization, the active enzyme centrifugation (AEC) 
[33] was performed. Briefly, MERS-CoV MP" of 15 pl (1 mg/ml) was added into the small well 
of the band-forming centerpiece before the cell assembled. Then 330 ul of peptide substrate at 
0, 200 and 400 uM in DO were respectively loaded into the bulk sample sector space. At a 
rotor speed of 42,000 rpm, the protein solution flowed into the substrate-containing channel 
and form a protein band. It can be detected by absorbance at 250 nm. During the centrifuga- 
tion, the sediment protein continuously met and cleaved the substrate, which can be detected 
by absorbance change at 405 nm. The dataset from the multiple scans at 250 nm at various 
time intervals were fitted to a continuous c(s) distribution model using the SEDFIT program 
[31], while the first five scans (0-30 min) at 405 nm were used to derive the product concentra- 
tion and then initial velocity values. 
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Isothermal titration calorimetry (ITC) 


The protocol followed that of previous studies [28] with some modifications. Apparent dissoci- 
ation constants and stoichiometry of the enzyme-ligand interactions were measured by a Ther- 
mal Activity Monitor 2277 from TA instruments (New Castle, DE). Calorimetric titrations of 
the peptide substrate TQ6-pNA (0.5 mM in a 250-ul syringe) and MP"® (6 uM in a 4-ml 
ampoule) were carried out at 25°C in 10 mM phosphate buffer (pH 7.6). The peptides were 
titrated into the enzyme using a 10-l aliquot for each injection with a time interval of 20 min. 
A control experiment in the absence of enzyme was performed in parallel to correct for the 
dilution of heat. The data obtained was then analyzed by integrating the heat effects normalized 
against the amount of injected protein using curve-fitting based on a 1:1 binding model. This 
involved the use of Digitam software (TA instruments, New Castle, DE). 


Results and Discussion 
Recombinant MERS-CoV MP” preparation 


As part of the present study, an expression vector was constructed and the BL21 (DE3) STAR 
(Invitrogen) strain of E. coli were used to express MERS-CoV MP"®. Unlike SARS-CoV MP"® 
[25, 28], the MERS-CoV MP” with 6 x His-tag retained at the C-terminus cannot be expressed. 
Instead, the bacteria are able to express the M?"° when there is a N-terminal 6 x His-tag fusion 
that can be removed during the purification. However, thrombin digestion leaves two extra res- 
idues (Gly-Ser) at the N-terminus of M?"°, resulting in protein with no proteolytic activity 
(data not shown). Therefore we used SARS-CoV papain-like protease [12, 30, 34], which is a 
highly active viral deubiquitinase and does not leave any residues at the N-terminus of MP". 
After gel-filtration, the purity of authentic N-terminus MP" was about 99% (S2 Fig). The size 
of the MERS-CoV MP”® was found to be close to 30 kDa, while any uncut protein was located 
at higher molecular weight position. The typical yield was about 10 mg after purification from 
0.8 liter of E. coli culture. 


Overall structure of MERS-CoV MP’? 


The structure of the MERS-CoV MP? was determined at 3.0 A resolution by X-ray crystallog- 
raphy (Table 1 and Fig 1A). The crystal packing belonged to space group C222,, with unit-cell 
parameters a = 87.2, b = 94.0, c = 155.1 Aanda= B=y = 90°. The final atomic model contain- 
ing two MP"° molecules in a crystallographic asymmetric unit agrees well with the crystallo- 
graphic data and the expected values of geometric parameters (Table 1). There are no residues 
in the disallowed region of the Ramachandran plot, while 81.3% of the residues are in the most 
favored region. 

The overall dimeric structure of MP"° is similar to that of SARS-CoV MP"; although a rela- 
tive shift of 10° to 30° could be observed for the two domain III within the dimer (S3A Fig). 
Indeed, the r.m.s. distance between equivalent Ca atoms of the domain I+I] of the two struc- 
tures is 0.9 A, while that between the domain III of the two structures is 3.1 A. Compared with 
the structures of the ligand-bound complex (PDB entry 4YLU) and C148A mutant (PDB entry 
AWME) [14, 15], the r. m. s. distance is 0.8 and 0.7 A over 540 Co atom pairs, respectively (S3B 
Fig). This indicates that the dimeric structures show no significant difference; although the 
present structure is a free enzyme and the other two structures involve enzyme-ligand com- 
plexes and higher resolution. Besides, there is minor difference between the present structure 
and bat-CoV HKU4 MP?" (PDB entry 2YNA), with 80% sequence identity, as the r. m. s. dis- 
tance over 539 Ca atom pairs is 0.8 A (S3B Fig). 
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Data Collection 


Dimerization and Catalysis of MERS-CoV Main Protease 


Space group C222, 
Cell dimensions 
a, b, c (A) 87.2, 94.0, 155.1 
a, B, y (°) 90, 90, 90 
Resolution? (A) 30-3.0 (3.11-3.0) 
Rinerge” (%) 17 (62.9) 
I/ol 10.2 (3.1) 
Completeness (%) 94.6 (95.2) 
Redundancy 5.5 (5.6) 
Refinement 
Number of reflections 11,675 (1,549) 
R factor® (%) 21.8 
Free R factor? (%) 28.2 
Number of atoms 
Protein 4,580 
Water 0) 
Average B-factors for protein atoms (A?) 81.6 
R.m.s. deviations 
Bond length (A) 0.009 
Bond angles (°) 1.5 
Ramachandran plot statistics (%) 
Most favored region 81.3 
Additional allowed region 13.6 
Generously allowed region 5.2 
Disallowed region 0 


* The numbers in parentheses are for the highest-resolution shell. 


ehaee = b> > hi = (I,) 
Ah i 


mean intensity of multiple corresponding symmetry-related reflections. 
SS 
h 


4 Free R is R calculated using a random 5% of data excluded from the refinement. 


iy Sole where |; is the integrated intensity of a given reflection and (/,) is the 
h i 


FF, oe F;, where Fp and Ff are the observed and calculated structure factors, respectively. 
h 


doi:10.1371/journal.pone.0144865.t001 


Interestingly, the dimerization interface situation with the MERS-CoV MP?" was found to be 
different to that of SARS-CoV MP"? where there are four amino acid pairs with intermolecular 
polar interactions (Serl.. .Glu166, Arg4. . .Glu290, Ser123. ..Arg298 and Ser139. . .GIn299). 
There are only two pairs of intermolecular hydrogen bonds, Ser1. . .Glu169 and Ser142. . .GIn299 
that are associated with the dimer surface of MERS-CoV M?” according to the current structure 
(Fig 1B). This led us to compare the dimerization and catalytic activity of the two types of MP". 
In the present study, in addition to using the wild-type MERS-CoV M?"°, we also mutated several 
residues at the dimerization interface in order to evaluate their role in dimerization and catalysis 
of MERS-CoV MP" (see below). 


Correlation between dimerization and catalysis of MERS-CoV MP'° 


To compare catalysis between the two M?"°, TQ6-pNA, a peptide substrate for SARS-CoV 
MP"® [25, 28], was used to measure the proteolytic activity. At first, the dependence of the initial 
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Fig 1. The structure of MERS-CoV MP". (A) The overall structure of the dimeric M?"°. The two protomers 
are shown as aribbon and as a surface model, respectively. The negatively and positively charged regions 
on the molecular surfaces are colored red and blue. The spheres indicate the catalytic dyad, His41-Cys148. 
(B) Stereo view of an overlay of the dimerization interface and the active site of MERS-CoV MP” (in cyan and 
orange) with that of SARS-CoV MP"? (grey; [16]). The red dashed lines indicate polar interactions between the 
two protomers of MERS-CoV M?"°, while the black dashed lines show polar interactions between those of 
SARS-CoV MP". The C atoms of the modeled substrate P4-P1 residues (from the structure of C148A mutant 
[15]) are colored magenta. The structural figures in this paper were produced using PyMol (http://www.pymol. 
org/). 


doi:10.1371/journal.pone.0144865.g001 


velocity on enzyme concentration was analyzed and showed a nonlinear upward correlation 
(Fig 2A). The pattern is similar to that of SARS-CoV MP", as the monomeric M?™ may not 
have catalytic activity [28]. However, MERS-CoV MP" displayed a sigmoid curve for its rate 
constant pattern at various substrate concentrations (Fig 2B, open circles); this contrast with 
SARS-CoV MP"®, which exhibited a classical saturation curve (S4 Fig). The results were then 
fitted to the Hill equation (Eq 2) in order to evaluate the kinetic parameters (‘Table 2). The kcat 
(2.33 s') of MERS-CoV MP” is close to that of SARS-CoV MP"? (2.11 s'), while the Hill con- 
stant was 1.8, suggesting a significant degree of positive cooperativity among the M?”® proto- 
mers. The comparable activity levels of the two MP" in the present study is dissimilar to the 
results obtained during a recent study in which the activity of the MERS-CoV MP?” was found 
to be 5-fold lower than that of SARS-CoV [14]. Using different substrates may cause the differ- 
ence. Tomar et al. [14] used a longer peptide substrate with residues present at both P and P’ 
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Fig 2. Proteolytic activity assay of MERS-CoV MP"? and its mutants. (A) Plot of initial velocities (v;) versus 
the concentration of wild-type MERS-CoV MP’°. The concentration of substrate (TQ6-pNA) was 600 uM. The 
line represented the best-fit results according to the nonlinear dependence equation [28]. (B) and (C) Plots of 
rate constant (kops) versus the concentration of peptide substrate for MERS-CoV MP"? (by circles), the V4R 
mutant (by lower triangles), the T126S mutant (by squares), the E169A (by diamonds), the M298R (by upper 
triangles) and the T126S/M298R (by hexagons) are indicated. The lines represented the best-fit results 
according to the Hill equation (Eq 2; wild-type, V4R, T126S and E169A mutants) or the Michaelis-Menten 
equation (Eq 1; M298R and T126S/M298R mutants). The protein concentrations of the wild-type, V4R, 
T126S, E169A, M298R and T126S/M298R mutants used for the assay were 0.3, 0.4, 0.7, 1.2, 0.15 and 

0.26 uM, respectively. All the assays were performed in 10 mM phosphate (pH7.6) at 30°C and repeated 
twice to ensure reproducibility and the error bars were shown. The kinetic parameters are shown in Table 2. 


doi:10.1371/journal.pone.0144865.g002 
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Table 2. The kinetic parameters and dissociation constants of MERS-CoV MP"°. 


Proteins Kinetic Parameters* Dissociation Constant? 
Km (uM) Keat (S') h No substrate(yM) With 600 uM substrate (uM) 
MERS-CoV Pr? i 
Wild-type - 2.33 + 0.13 1.8 + 0.04 7.7+0.3 0.7 + 0.04 
V4R - 0.96 + 0.05 2.7+0.2 23.0 + 0.4 15.2+0.3 
T126S - 0.56 + 0.04 2.0+0.2 33.7 + 0.9 13.9 + 0.1 
E169A - 0.41 + 0.02 2.140.1 14.3+0.2 14.1+0.5 
M298R 181.0 + 24.0 7.91 0.49 - 1.1401 0.7 +0.01 
T126S/M298R 419.4 + 63.9 4.63 + 0.37 - 2.8 + 0.1 0.9+0.01 
SARS-CoV Mpro 890 + 130 2.11 40.15 - 0.7 + 0.02 1.7 + 0.03° 


@ Kinetic data of SARS-CoV MP"? and MERS-CoV M298R and T126S/M298R mutants were fitted to the Michaelis-Menten equation (Eq 1), while those of 
the others were fitted to the Hill equation (Eq 2). The Reg, were from 0.985 to 0.999, respectively. All the assays were repeated twice and the average 


values were used for the fitting. 


> The values were derived from a global fit of the AUC data to a monomer-dimer self-association model by SEDPHAT [32]. The experiments for the assay 
were obtained at protein concentration of 1.5 to 30 uM. 


° The value was from our previous studies for comparison [28]. 


doi:10.1371/journal.pone.0144865.t002 


site. However, in the present studies we used a peptide substrate that contains only P site resi- 
dues. Besides, they utilized FRET substrate and could only be used at low substrate concentra- 
tions to prevent the inner-filter effect, while we are able to use higher substrate concentrations 
to capture the kinetic parameters. 

The cooperativity phenomenon associated with the MERS-CoV MP" is similar to that of 
the SARS-CoV MP" R298A/L monomer mutants; these were found to show monomer to 
dimer conversion during catalysis [28]. As a result of the above, we investigated the quaternary 
structure of the MP"° by AUC (Fig 3). The cumulative spectra (Fig 3A) were analyzed using the 
continuous c(s) distribution model and the results suggested that MERS-CoV MP"° is a mono- 
mer in phosphate buffer (S5 Fig) and this contrasts with a distribution of 30% monomer and 
70% dimer in the presence of 600 uM TQ6-pNA (Fig 3B). We also measured the size distribu- 
tion of MP" at various TQ6-pNA concentrations (S6 Fig). The results indicated that the sedi- 
mentation coefficient of the major species was shifted as the substrate dosage changed (SOB 
Fig). More substrate led to the major species moving close to the dimer position. 

However, before the centrifugation, the enzyme had been mixed with the substrate and the 
catalysis began, resulting in a mixture of substrate and product with enzyme. It is unable to 
confirm that our observation is a substrate-induced or substrate/product-induced dimeriza- 
tion. To solve this, a modified AUC technique, AEC [33], was utilized to detect the quaternary 
structure change in the absence and presence of substrate (Fig 4). Here the enzyme solution 
was put into the small well of band-forming centerpiece and then flowed into the substrate- 
containing channel when the centrifugation began. During the centrifugation, the protein layer 
gradually sediment and continuously met peptide substrates. Not surprisingly, there was a 
broad distribution between the monomer and dimer species in the presence of 200 uM peptide 
substrate, while a major species shifted to the dimer in 400 uM substrate (Fig 4B). It suggests 
that MERS-CoV MP?" acts as a rapid self-associated and substrate-induced dimerization. Using 
different strategies, Tomar et al. [14] confirmed that inhibitor binding can also induce and 
maintain the dimerization of M?"®. On the other hand, we measured the velocity of the product 
formation during the centrifugation; although the rate (0.017 M/s) is 10-fold lower than that 
by the spectrometric assay (Fig 4B). 
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Fig 3. The continuous size distribution change of MERS-CoV MP"°. (A) Typical trace of absorbance at 
250 nm of the enzyme during the SV experiment. The protein concentration was 0.5 mg/ml. For clarity, only 
every fifth scan is shown. The symbols represent the experimental data and the lines represent the results 
fitted to the Lamm equation using the SEDFIT program [31]. (B-G) Continuous c(s) distribution of wild-type, 
V4R, T126S, E169A, M298R and T126S/M298R mutants. The distributions in D20 containing 10 mM 
phosphate (pH 7.6) are shown by solid lines and those in the same buffer but with 600 UM TQ6-pNA 
substrate are shown by dashed lines. The left vertical dotted line indicates the monomer position and the right 
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dotted line represents the dimer position. The residual bitmap of the raw data and the best-fit results are 
shown in the insets. 


doi:10.1371/journal.pone.0144865.g003 


To quantitatively characterize the monomer-dimer equilibrium of M?"® in the absence and 
presence of substrate, the AUC results at protein concentrations of 1.5, 6 and 15 uM were glob- 
ally fitted to the monomer-dimer self-association model (Table 2). The analysis indicated that 
the Kg value of MERS-CoV MP" in the absence of substrate was 7.7 uM and this decreased 
11-fold in the presence of substrate, which brings it close to the value for SARS-CoV MP"° 
(1.7 uM; Table 2). Thus it can be concluded that, like the SARS-CoV MP*° R298A mutant [25], 
the presence of substrate is able to induce the dimerization of MERS-CoV MP". In addition, 
although the Kg values of wild-type SARS-CoV MP?"® without or with substrates show no signif- 
icant difference (Table 2), it was possible to detect substrate-induced dimerization at a protein 
concentration of 1 uM by AEC [33]. 

Previous studies have demonstrated that a conserved residue, Glu166, plays a pivotal role 
in connecting the substrate binding site with the dimerization interface of SARS-CoV MP*° 
[28]. Here the equivalent residue, Glu169, was mutated and its influence on MERS-CoV MP° 
evaluated (Fig 1B). Not surprisingly, compared to the wild-type enzyme, the activity (k,at) of 
E169A showed a 6-fold decrease (Fig 2B, open diamonds and Table 2). Furthermore, AUC 
analysis suggested that this enzyme consisted of a single species close to monomeric form 
even in the presence of substrate (Fig 3E). The Kg values of the E169A mutant without sub- 
strate was 14 uM, which is 2-fold higher than that of the wild-type and this did not decrease in 
the presence of substrate (Table 2). The results suggest that mutation of Glu169 is able to 
block substrate-induced dimerization and that this results in a decrease in enzyme activity. 
Unexpectedly, as well as the monomer, some octamer (6.3%) with a sedimentation coefficient 
of 4.9 S was observed in the presence of substrate (Fig 3E). Previous studies have suggested 
that a super-active octamer of SARS-CoV MP” can be locked by 3D domain swapping [35]. 
In this study, the presence of the octamer form of the MERS-CoV MP" E169A mutant in the 
presence of substrate may explain why this mutant is not totally inactive. Taking the above as 
a whole, the dual role of the conserved Glu residue in catalysis and dimerization is consistent 
for both MP". 

Our results confirmed that there are fewer intermolecular polar interactions at the dimeriza- 
tion interface of MERS-CoV MP” and this results in the enzyme being in the monomer form 
in aqueous buffer; this contrasts with SARS-CoV MP", which is mostly in the dimer form in 
similar circumstances due to the greater number of intermolecular interactions (Figs 1B and 
3B). Based on the sequence alignment, three residues in the dimerization interface vary in coro- 
naviral MP"® sequences (S1 Fig). In SARS-CoV MP", residues Arg4, Ser123 and Arg298 are 
different to the equivalent ones in MERS-CoV MP", which are Val4, Thr126 and Met298, 
respectively. Based on the above, three single mutants, V4R, T126S and M298R, were generated 
and their proteolytic activity and dimerization were assessed. Unexpectedly, both the V4R and 
the T126S mutant showed a higher Kg for the dimer to monomer and a lower level of activity 
than wild-type MP"° (Fig 2B and Table 2). Kg values were decreased by 1.5-fold to 2.4-fold for 
the two mutants in the presence of substrate, which suggests reduced substrate-induced dimer- 
ization (Fig 3C and 3D). Based on the current structure, residue Val4 showed a hydrophobic 
contact with another protomer’s residue Gly141 (Fig 1B). Mutation of valine to arginine may 
lose the contact. Furthermore, the side chain of residue Val131 of MERS-CoV MP” is hydro- 
phobic while that of the equivalent residue, Cys128 of SARS-CoV MP", is hydrophilic (Fig 1B). 
Val131 is close to the Arg4 and will disfavor the electrostatic interaction of Arg4. . .Glu290 in 
MERS-CoV M?"’. These variance may result in V4R failed to form a stable dimer. For the 
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Fig 4. Charactering the substrate-induced dimerization of MERS-CoV MP"? by AEC [33]. (A) A trace of 
absorbance at 250 nm for the enzyme in the presence of 200 uM of TQ6-pNA during the experiment. For 
clarity, only every twice scan and every fifth measuring points per scan are shown. The protein of 15 yg in 10 
mM phosphate (pH7.6) was used. The substrate was dissolved in D2O to give a higher density, which can 
maintain the protein band when the centrifugation begins. The symbols are experimental data and the lines 
are the results fitted to the Lamn equation using SEDFIT program [31]. The continuous size distribution of the 
best-fit result was shown in panel B colored by green. (B) Continuous c(s) distributions of MERS-CoV MP” at 
the concentration of TQ6-pNA of 0, 200 and 400 uM. The labels M and D showed the position of the monomer 
and dimer species. (C) Monitoring the enzyme activity of MERS-CoV MP"? during the centrifugation. The 
same cell in panel A was followed and the absorbance at 405 nm trace for the released product (pNA) was 
shown. We used the spectra of the first 30 min to calculate the product concentration at different time. The 
inset plot showed the product versus time and the slope of the line represented the initial velocity. 


doi:10.1371/journal.pone.0144865.g004 
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T126S mutant, after compared with the other two MERS-CoV MP” structures (PDB entry 
4YLU and 4WMBE), we found that the side chain of Thr126 is free of rotation. It is able to make 
a hydrophobic interaction with the residue Tyr121, resulting in a 180°-rotation of the phenol 
ring, and lead to a hydrogen-bond with the backbone amide of Leu144 (Fig 1B). It further 
results in the side chain of Leul44 toward another protomer’s Ile300 and make a hydrophobic 
contact for the two protomers. Mutation of Thr126 to serine will lose this contact and may dis- 
favor the dimerization. 

By way of contrast, the M298R mutant resulted in a stable dimer form in phosphate buffer 
(Fig 3F). The Kg value for the mutant in the absence and presence of substrate were 1.1 and 
0.7 4M, respectively, which are very close to those for SARS-CoV M?"® (Table 2). Furthermore, 
the stable dimer form showed higher proteolytic activity and the mutation transformed the 
enzymes rate constant pattern at various substrate concentrations into a classical saturation 
curve (Fig 2C, open triangles). In addition, the K,, of the M298R mutant is 5.8-fold lower than 
that of SARS-CoV MP", which suggests a higher substrate binding affinity (Table 2). It can be 
concluded that mutation of residue Met298 to arginine within the MERS-CoV MP” results in 
the stabilization of the dimer formation, which in turn gives rise to more efficient catalysis. 
Based on our structure, Arg298 is able to make a hydrogen bonding interaction with Thr126. 
Residue Thr126 can be replaced by serine because the T126S/M298R double mutant also 
shows similar dimerization characteristics (Fig 3G) and saturated catalytic pattern to that of 
the M298R mutant (Fig 2C, open hexagons and Table 2). However, it can only be achieved in 
the presence of Arg298, not Met298. 


Substrate binding analysis of MERS-CoV MP” by ITC 


The sigmoid nature of the curve describing the rate constant pattern at various substrate con- 
centrations (Fig 2B) means that it is not possible to obtain a K,, value for this enzyme, which 
would allow us to evaluate the substrate binding affinity of the wild-type MERS-CoV MP" and 
its mutants; the exceptions being the M298R single mutant and the T126S/M298R double 
mutant (Fig 2). To further delineate the binding of substrate to the enzyme, ITC was used to 
measure the Kg for the substrate (or substrate/product)-enzyme complex and the binding stoi- 
chiometry (N) (Fig 5). During the titration, the enzymatic hydrolysis might produce additional 
heat, resulting in higher AH. So we can only compared the N and Kg of the wild-type M?"° with 
those of E169A and M298R mutants. The three enzymes exhibited similar N (0.89 to 1.06) and 
Kg (14.2 to 20.3 uM). This suggested that the monomeric and dimeric MP" show quite the 
same substrate binding affinity. Such phenomenon is also found in monomeric and dimeric 
SARS-CoV MP"®, whose N and Kg for the same substrate were 0.97 to 1.05 and 29.9 to 33.8 uM, 
respectively [28]. With the Ky, ((k.y+keat)/ki), Kear and Kg (k_1/k,), we are able to calculate the 
k, and k ,, the rate of the association and dissociation of the enzyme-substrate complex. The 

k, and k., for the M298R mutant and substrate is 0.049 s'uM” and 0.98 s’, while those for 
SARS-CoV MP" is 0.00246 s'M”' and 0.084 s’. The rate constants for M298R mutant 
showed 20- and 12-fold higher than those for SARS-CoV MP"°. The more rapid association 
and relatively slower dissociation of enzyme-substrate complex may be used to explain why the 
catalytic efficiency (k.at/Km) of the MERS-CoV MP" M298R mutant is higher than that of 
SARS-CoV MP" (Table 2). 


Conclusion 


The crystal structure of authentic N-terminus MERS-CoV MP" was determined and this was 
found to involve a dimeric form with less intermolecular polar interactions. Biochemical and 
AUC studies indicated that MERS-CoV M?"® shows almost the same proteolytic activity as 
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Fig 5. ITC titration of peptide substrate binding to MERS-CoV MP’? wild-type (A), E169A (B) and M298R 
(C) mutants. Raw data are shown in the top panel of each panel and represent the power input into the 
sample ampoule over time. Integrated data are shown in the bottom panels in terms of the total energy 
required for equilibration as a function of the molar ratio of substrate to M?°. The solid circles show the 
observed values and the lines represent the fitted results by ligand binding analysis using the Digitam 
program (TA instruments). The best-fit parameters for stoichiometry (N) and Kg are shown in the bottom 
panels. 


doi:10.1371/journal.pone.0144865.g005 


SARS-CoV MP"; although it is a monomer in aqueous buffer and displays substrate-induced 
dimerization (Fig 6). A conserved residue, Glu169, plays an essential role in the substrate- 
induced dimerization of both MERS-CoV MP"° and SARS-CoV MP". Moreover, mutation of a 
residue in the dimerization interface, M298R, was found to result in a more stable dimer form 
in aqueous buffer that had higher enzyme activity; while other two mutations, V4R and T126S, 
showed the reverse effect. Critical assessment of the residues important to dimerization of and 
catalysis by MERS-CoV MP"® provides valuable insights into the mechanism that controls the 
monomer-dimer switch of important and valuable enzyme. 
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Fig 6. Schematic model of the catalytic process of MERS-CoV MP"°. Substrate binding results in an 
active protomer (square) and then dimer formation. After a catalytic cycle, the dimer will dissociate to form the 
inactive protomers (circle). The M298R mutant, like SARS-CoV MP", is able to maintain a dimer form without 
dissociation before the next round’s substrate binding and catalysis. 


doi:10.1371/journal.pone.0144865.g006 
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Supporting Information 


S1 Fig. Sequence alignment of betacoronaviral M?*°. Modified from an output from ESPript 
[36]. The green ovals indicate the catalytic dyad, while the blue ovals indicate the residues mak- 
ing intermolecular polar contact in the dimer interface of SARS-CoV Mpro. Magenta oval 
indicate the residue Glu playing dual role for the dimer interface and substrate binding site. 
Accession numbers are as follows: MERS-CoV, NC_019843.2; Bat-CoV_HKU4, ABN_010865.1; 
HCoV_HKU1, NC_006577.2; MHV_A59, NP_068668.2; SARS-CoV, NP_828863.1. 

(PDF) 


$2 Fig. Expression and purification of recombinant MERS-CoV MP”. Protein identification 


by SDS-PAGE. M: molecular marker. Lane 1-5: cytoplasmic fraction, flow-through, elute from 
the nickel affinity column, flow-through after 4-h’s PL?” treatment and protein fraction from 
S-300 gel-filtration column. 

(PDF) 


$3 Fig. Comparison with the structure of MERS-CoV MP” with that of SARS-CoV MP*® 
(A) and those of ligand-bound complex, dimeric C148A mutant and bat-CoV HKU4 MP”° 
(B). (A) An overlay of the current structure of MERS-CoV MP" (cyan and orange) with that of 
SARS-CoV (grey; PDB entry 1uk3). The red arrows show the orientation change affecting the 
two domain IIIs. Spheres show the two catalytic dyads. (B) Overlay of the current structure 
with ligand-bound complex (red; PDB entry 4YLU), dimeric C148A mutant (magenta; PDB 
entry 4WME) and bat-CoV HKU4 MP” (yellow; PDB entry 2YNA). 

(PDF) 


S4 Fig. Activity assay of SARS-CoV M?"®. The plot of rate constant (k,,,) versus the concen- 
tration of TQ6-pNA are indicated. The line represented the best-fit results according to the 
Michaelis-Menten equation (Eq 1). The protein concentration was 1.1 uM. The assays were 
performed in 10 mM phosphate (pH7.6) and repeated twice to ensure reproducibility and the 
error bars were shown. The kinetic parameters are shown in Table 2. 

(PDE) 


$5 Fig. Continuous size distribution of MERS-CoV MP” at various protein concentrations. 
The protein concentration of 1.5 (solid line), 6 (dotted line) and 15 uM (dashed line) of 
MERS-CoV MP"° were used and monitored the size distribution by AUC. The best-fit results 
suggest that the major species was a monomer with minor shift of sedimentation coefficient 
(1.74 to 1.97), while the calculated molar mass was from 35.8 to 36.6 kDa. The residual bitmaps 
of the raw data and the best-fit results are shown in the insets. 

(PDF) 


S6 Fig. Effect of substrate concentration on the dimerization of MERS-CoV M?”. (A) Con- 
tinuous c(s) distribution of the enzyme at peptidyl substrate (TQ6-pNA) concentrations of 

0 uM (solid circles), 80 uM (open circles), 200 uM (solid triangles), 340 uM (open triangles), 
400 uM (solid squares) and 600 uM (open squares). The protein concentration was 0.25 mg/ 
ml. (B) Sedimentation coefficient shifts of the major species of MERS-CoV MP?” at different 
TQ6-pNA concentrations. 

(PDF) 
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